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Abstract 

We have studied the decay B — > Div, where i = e or \x. From a fit 
to the differential decay rate dT/dw we measure the rate normalization 
•Fd(1)|V^| and form factor slope p 2 D) and, using measured values of t b , 
find T(B Div) = (12.0 ± 0.9 ± 2.1) ns -1 . The resulting branching 
fractions are B(B° -> D+£~u) = (1.87 ± 0.15 ± 0.32)% and B(B~ -> 
D°£~i>) = (1.94 ± 0.15 ± 0.34)%. The form factor parameters are in 
agreement with those measured in B — > D*iv decays, as predicted by 
heavy quark effective theory. 
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Exclusive semileptonic B meson decays provide information about both the weak and 
strong interactions of quarks. The rate for these decays is proportional to the square of the 
CKM matrix element \V c b\ JT|, while the dynamics of these decays, as expressed in the decay 
form factors, provide information about the QCD potential which binds quarks together 
as hadrons Heavy quark effective theory (HQET) J3J has made it possible to extract 

values of \V c b\, with relatively little model dependence, through a measurement of the decay 
rate at the point of zero recoil of the daughter meson. The decay B — ► D*lv has provided 
measurements of the decay rate at the point of zero recoil and the form factor slope, giving 
very accurate values of \V c b\ [|5| and information about the shape of the Isgur-Wise function 
that describes the dynamics of heavy-to-heavy meson transitions. 

The pseudoscalar-to-pseudoscalar decay B — > Dtv can provide the same information, 
although it will be less precise because of the smaller overall decay rate for this mode than 
for B — > D*tv, the smaller rate near the point of zero recoil, and 0(1/ Mq) corrections to 
the decay rate at that point that are not present in B — > D*tv. However, this mode should 
yield a value of \V c i\ that is consistent with other measurements, and HQET predicts 
that the form factor parameters should be very nearly the same as in the B — > D*tv decay. 
In this Letter, we present a measurement of the differential decay rate dT/dw in the decay 
B — > D£u, from which we extract the decay rate normalization at the point of zero recoil 
and the form factor slope. The variable w = (Mg + M|> — g 2 )/(2MeMo), where q 2 is the 
invariant mass squared of the lepton-neutrino system, is the kinematic variable of HQET, 
and is equal to the relativistic 7 factor for the D meson in the B meson rest frame. From 
these results and the measured B lifetime, we obtain the partial width for the decay and 
convert it to branching fractions. 

This study is based on an T(AS) data sample of 3.16 fb _1 (3.34 xlO 6 BB pairs) accumu- 
lated by the CLEO experiment at the Cornell Electron Storage Ring (CESR). The CLEO 
detector contains three concentric wire chambers that detect charged particles and a 
CsI(Tl) electromagnetic calorimeter that detects photons, all within a 1.5 T superconduct- 
ing solenoid. 

The undetected neutrino complicates analysis of semileptonic decays. Using the her- 
meticity of the CLEO detector, we reconstruct the neutrino by inferring its four-momentum 
from the missing energy (E miss = 2£ beam - £i?i) and missing momentum (p miss = ~T,Pi) 
in each event, where Ei and pi are the energy and momentum of each detected particle i in 
the event, as was done in the CLEO measurement of exclusive b — > ulv decay rates ||. In 
the process e + e _ — > T(AS) — > BB, the total energy of the beams is imparted to the BB 
system. At CESR, that system is at rest, so the neutrino combined with the signal lepton 
and D meson should satisfy the energy constraint AE = (E p + Eg + E D ) — E heiim = and 
the momentum constraint M cand = [-Ebeam — \p„ + Pi + Pd| 2 ] 5 = M B . We select candidates 
with 5.2650 < M cand < 5.2875 GeV and -100 < AE < 500 MeV; the requirement on AE is 
asymmetric about zero to reject feeddown from B — > D*lv decays, which have AE values 
of about -150 MeV when reconstructed as B — > Dtv decays. 

To suppress events in which p miss misrepresents p p , we reject those events with multiple 
leptons or a total charge more than one unit from zero because they indicate other missing 
particles. We further require that M^ iss = -E^ iss — |p m iss| 2 f° r eac h event be consistent 
with zero. Surviving signal events show a resolution in |p m i SS | of 110 MeV/c. Because the 
resolution on -E m i SS is about 2.1 times larger, we take (E p ,p p ) = (|p m iss|;Pmiss)- 
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Information from calorimeter and tracking measurements including specific ionization is 
combined to identify electrons with p > 600 MeV/c over 90% of the solid angle. Particles 
are considered muons if they register hits in counters deeper than 5 interaction lengths over 
the polar angle range | cos#| < 0.85. Candidate leptons must have 0.8 < pi < 2.4 GeV/c, 
where the lepton identification efficiency averages more than 90%; the probability that a 
hadron is misidentified as an electron (muon), a "fake lepton", is about 0.1% (1%). 

The leptons and neutrinos are then combined with D mesons, which are identified in 
the decay modes D° —>■ K~n + and D + — ► K~tt + tt + . (The charge conjugate mode is always 
implied.) Hadron mass assignments are made by requiring that the kaon and lepton have 
the same charge. To reduce large backgrounds to D + — ► K^-k + -k + decays from random 
track combinations, we require that the D + daughter pions (kaon) have measured specific 
ionization values consistent with the assumed particle hypothesis within 3.5 (3) standard 
deviations. D° candidates are required to have invariant mass satisfying 1.850 < M^n < 
1.880 GeV, and D+ candidates to satisfy 1.855 < M Kn7r < 1.885 GeV. This ±15 MeV range 
is about twice the experimental resolution on the invariant mass. Since the D must come 
from a B decay, we require po < 2.6 GeV/c. We also require that the angle between the 
directions of the D candidate and the lepton-neutrino system have a physical value;the angle 
is calculated from the D meson and lepton momenta and the beam energy only, without 
examining the missing momentum. 

Backgrounds arise from e + e~ — > qq/r + r~ (continuum), fake leptons, random track com- 
binations that form D candidates, feeddown from other B — > DXiv decays, and random 
combinations of D candidates and leptons from different parent B mesons. The continuum 
backgrounds are reduced by requiring that the ratio of Fox- Wolfram moments H2/H0 
be less than 0.4. Backgrounds from B — > D*X£v decays are reduced by eliminating events 
that include Dn or D7 pairs that are consistent with D* decay. We find 303 D°£~v and 
714 D + £~v candidates that satisfy all of these requirements. The average reconstruction 
efficiency, as determined by a Monte Carlo simulation of the CLEO detector, is 2.82% for 
B~ — > D°£~u events and 2.57% for B° — > D + £~u events. Figure [I] shows the M can d distri- 
bution observed in the data. 

We estimate the continuum background using data collected 60 MeV below the T(4S) 
energy and the fake lepton background by applying measured fake rates to nonleptonic data. 
We estimate backgrounds from random track combinations that form D candidates using 
events in sideband regions on either side of the K~tt + (tt + ) invariant mass signal region under 
the assumption that the magnitude of this background is linear in K~tt + (tt + ) invariant mass. 
This is the largest background in the D + £~v sample. 

Feeddown backgrounds are modeled through Monte Carlo simulations, using an event 
generator that accounts for all angular correlations among the decay products, and a full sim- 
ulation of the CLEO detector. The magnitude of the B — > D*£u background is normalized 
to the measured rate for this decay ]I0| . This is the largest background in the D°£~u sample, 
as the kinematics of the signal and background decays are so similar, and because D* mesons 
are more likely to decay to a D° than to a D + . B — ► D**£u processes, where D** represents 
a variety of charm mesons with radial and angular excitations and nonresonant D^tc states, 
are modeled with the ISGW2 and Goity and Roberts |IT| models. We normalize these 
backgrounds to a set of rates Jl2]| that are consistent with existing measurements [13| and 
the total B semileptonic decay rate and lepton momentum spectrum [F4]. 
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FIG. 1. M can( j distributions for the D°£~i> (top) and D + l~v (bottom) modes. The points are 
the data, the shaded component is the continuum and fake lepton background, the diagonal hatch 
is the combinatoric background, the vertical hatch is the B — > D*lv feeddown, and the Crosshatch 
is the B — > D**lv and other backgrounds. The unshaded area is the prediction of a Monte Carlo 
simulation of the signal, normalized to the measured decay rate. 
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TABLE I. Event yields and background estimates. Errors are statistical only. 





D°(.-i> 


D+tTv 


Total Yield 


303.0 ± 17.4 


714.0 ± 26.7 


Continuum 


17.3 ± 5.8 


46.2 ± 9.4 


Combinatoric 


26.1 ± 4.6 


256 5 ± 10 8 

— ' J \.J • \.J — 1 — J- V 7 ■ V.- 1 


Fake lepton 


2.9 


2.3 


B -> D*£u 


107.2 


62.6 


B -» D**£u 


5.3 


9.9 


Other 


-0.1 


13.6 


Corrected Yield 


144.3 ± 18.9 


322.8 ± 30.3 



Remaining backgrounds, such as those from random combinations of D mesons and 
leptons and those from misreconstructed D mesons not accounted for in the combinatoric 
background estimate, are modeled by a Monte Carlo simulation that reproduces the general 
features of 5-meson decay, including the inclusive-lepton and D momentum distributions. 
The magnitudes of the various backgrounds are summarized in Table |. Figure |1| shows 
the M can d distributions of the backgrounds. After accounting for these backgrounds in our 
sample, the lepton momentum and decay angle distributions for selected events are consistent 
with B — ► Dlv decays. 

To extract the form factor parameters and partial width, we construct the distribution 



dT 



dw AN- 



T(45) 



N (w) + N+(w) 



r B -B D oe (w) T S oB D +e + (w) 



where A^x(45) is the number of T(45) events, e(w) is the reconstruction efficiency as a function 



of w, tb is the B lifetime, and B is the appropriate D — > K~7r + (-n- + ) branching fraction [10 
By combining the two decay modes in this fashion and assuming that the partial widths 
of the two modes are equal, the results are independent of the division of T(AS) decays 
between B + B~ and B°B°. This distribution is shown in Figure |2|a. 
The dT /dw distribution for B — > Div decays is || 

dT G 2 F \V cb \ 2 



., y m B + m D ) 2 {m D \/w 2 - lfF D {w) 2 . 
dw Attn 6 

The function Td{w) is the decay form factor, which can be parameterized by the expression 

Fd(w) = ^d(1)(1 - p 2 D (w - 1) + c D (w - I f). 

In the limit of infinitely heavy quarks, J-d{w) becomes the Isgur-Wise function £(w), and 
HQET predicts that J-"d(1) — 1; for finite mass quarks, ^0(1) can be estimated in the 
framework of HQET |Q. The values of p 2 D and Co are unknown; many previous form factor 
measurements have set cd = 0. We follow this convention, but explore our sensitivity to this 
assumption using various models that provide relations between p 2 D and cd- We perform 
a x 2 fit °f our dT /dw distribution to the convolution of the expected form with a function 
that accounts for detector resolution in w (±0.015). We allow two free parameters, the 
rate normalization •/ c d(1)H / c&| an d the form factor slope p 2 D . By integrating the fit over the 
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FIG. 2. (a) dT/dw distribution from the data (points), and the result of the c = fit to the 
distribution (histogram), (b) Measured values of .Fd(iu)|Vc&| (points), the result of the fit (solid 
line) along with its statistical errors (dashed lines), and the function Tb* (w)\ V c b\ obtained from 
an analysis of B — ► D*lv decays (dotted line). 
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TABLE II. Contributions to the systematic error (%) in the fit parameters and partial width. 
Simulation of the detector and the second B meson contribute to v simulation. 



Source 


\Vcb\ 


Pi 


r 


v simulation 


11.4 


12.5 


14.8 


Background normalization 


7.6 


15.4 


5.4 


B — > D*£v form factors 


1.8 


7.1 


2.0 


tb 


2.0 


2.5 


3.8 


Lepton ID 


1.0 




2.0 


K/ir ID 


0.1 


0.7 


0.8 


Luminosity 


0.9 




1.8 


D branching fractions 


2.8 




5.6 


Total 


14.3 


21.2 


17.5 



TABLE III. Results for various parameterizations of J-£>{w). The first errors are statistical 
and the second are systematic. p D and cd are entirely correlated in the Caprini-Neubert [17| and 
Boyd EH models, and fixed in the ISGW2 model §]. 



Model T D (1)| V cb |/10- 



Pl 



CD 



T (ns" 



X 2 /dof 



c D = 3.37 ± 0.44 ± 0.48 0.59 ± 0.22 ± 0.12 12.0 ± 0.9 ± 2.1 11.5/10 

c D free 4.57 ± 1.10 ± 0.92 1.84 ± 0.81 ± 0.53 1.75 ± 1.15 ± 0.56 12.3 ± 1.0 ± 2.2 10.3/9 

C-N 3.90 ± 0.65 ± 0.68 1.18 ± 0.37 ± 0.23 0.78 ± 0.27 ± 0.17 12.1 ± 1.0 ± 2.2 10.8/10 

Boyd 3.71 ± 0.60 ± 0.61 1.05 ± 0.38 ± 0.22 0.94 ± 0.43 ± 0.25 11.9 ± 0.9 ± 2.2 11.1/10 

ISGW2 3.25 ± 0.13 ± 0.27 0.64 0.61 12.1 ± 1.0 ± 2.0 11.6/11 



entire range of w, we obtain the partial width for the decay. In the fit, shown in Figure |2|a, 
X 2 = 11.5 for 10 degrees of freedom, and the correlation between the two parameters is 0.95. 

Systematic errors, summarized in Table [HJ are dominated by uncertainty in the decay 
model of the non-signal B and inaccuracies in detector simulation. These effects are in- 
vestigated by varying the fraction, charm semileptonic decay rate, charged-particle and 
photon-finding efficiencies, false charged-particle and photon simulation, charged-particle 
momentum resolution, and photon-energy resolution ||15|| . Uncertainties in the feeddown 



background normalizations and the B — > D*lv decay form factors [[U| have their most sig- 
nificant effect on the form-factor slope. Table [ITT] gives results for various models of Td{w)- 
The partial width is not sensitive to the choice of form-factor parameterization, but the val- 
ues of .Fd(1)|V^|, p D , and &d are sensitive to this choice. We use these results to determine 
our model uncertainties. 
Our final results are 

f D (l)\V cb \ = (3.37 ±0.44 ±0.48+°;^) x 10~ 2 
p 2 D = 0.59 ± 0.22 ±0.12+° 59 
T(B -> Dlv) = (12.0 ± 0.9 ± 2.1) ns^ 1 , 

where the first two errors are statistical and systematic, and the third arises from the form- 
factor model variations. This partial width leads to branching fractions of 

B(B° -> D+fv) = (1.87 ± 0.15 ± 0.32)% 
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B{B~ -> D°£-j?) = (1.94 ± 0.15 ± 0.34)%, 

where the errors are completely correlated between the two branching fractions. We obtain 
consistent results when the two decay modes are treated separately. Taking Fd(X) = 0.98 ± 
0.07 13, we find \V cb \ = (3.44 ± 0.45 ± 0.49lg;^ ± 0.25) x 10" 2 , where the last error arises 



from the uncertainty in .Fd(I). This value of \V c b\ is consistent with those obtained by other 
means |5|JT9|1. Figure ||b shows the measured values of J r D{w)\V c b\ and the result of the fit, 
along with the function -^d*^)!^! as measured in B — > D*tv decays by CLEO 0. In 
comparing these two measurements, we find JF D (1)/JF D »(1) = 0.96 ± 0.20 and p 2 D — p 2 D * = 
—0.25 ± 0.29. As predicted by HQET, the two form factors have similar normalizations and 
slopes. 

In summary, we have measured the rate normalization and form-factor slope in B — > Div 
decays. The resulting partial width leads to branching fractions for the charged and neutral 
B decay modes. The measured form factor parameters are consistent with those measured 
in B — > D*lv decays, as predicted by HQET, and the value of \V c b\ is consistent with other 
measurements. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Heisenberg Foundation, the Alexander von 
Humboldt Stiftung, Research Corporation, the Natural Sciences and Engineering Research 
Council of Canada, and the A. P. Sloan Foundation. 
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